Observation of Andreev Reflection Enhanced Shot Noise 
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We have experimentally investigated the quasiparticle shot 
noise in NbN/MgO/NbN superconductor - insulator - super- 
conductor tunnel junctions. The observed shot noise is signif- 
icantly larger than theoretically expected. We attribute this 
to the occurrence of multiple Andreev reflection processes in 
pinholes present in the MgO barrier. This mechanism causes 
the current to flow in large charge quanta {Andreev clusters), 
with a voltage dependent average value of m ~ 1 -I- |^ times 
the electron charge. Because of this charge enhancement ef- 
fect, the shot noise is increased by the factor m. 

PACS numbers: 72.70.-fm, 74.50.+r, 74.80.Fp 



A dc current / flowing through a vacuum tube or a tun- 
nel junction generates shot noise, time dependent fluctua- 
tions of the current due to the discreteness of charge carri- 
ers. Shot noise studies provide information on the nature 
of conduction not obtainable by conductance studies, e.g. 
the electric charge of carriers or the degree of correlation 
in the conducting system. For an uncorrelated system in 
which the electrons do not interact, the passage of car- 
riers can be described by a Poisson distribution. The 
spectral density of current fluctuations Si then equals 
full shot noise: S,j — 2ql = Ppoiaaon for zero frequency 
and temperaturetl. The charge quantum q is normally 
the electron charge e. 

In superconductor - normal metal (SN) systems An- 
dreev reflection occurs, causing an effective charge to be 
transferred of 2e. Due to this doubling of the charge, the 
shot noise in such a system is Modicted to have a maxi- 
mum of twice the Poisson noiseoQ. More recently, giant 
shot noise in the supercurrent is predicted in a single- 
channel superconductor - normal metal - superconductor 
(SNS) point contaciD which is attributed to transport of 
large charge quanta (q ^ e) at finite voltages caused by 
Multiple Andreev Refiection (MAR)l3. Observation of en- 
hanced charge quanta in SN or SNS structures requires 
a combination of conductance and shot noise measure- 
ments. Despite extensive theoretical work, experimental 
results in this field are rare. A recent experimentQ is per- 
formed on a NbN/c/Nb structure in which c is assumed 
to be a Nb constriction with a length of 7 nm and a di- 
ameter of 15 nm. At 9.5 K the structure acts like an 
NS interface but doubled shot noise is not observed. The 
predicted giant supercurrent shot noise is not observed 
either (at 4.2 K). 



From an applied point of view, shot noise in super- 
conducting structures is of interest since this noise forms 
a major limitation to the sensitivity of NbN/MgO/NbN 
Superconductor-Ipsulator-Superconductor (SIS) THz ra- 
diation detectorga. 

For these reasons, we have investigated quasiparticle 
current transport and shot noise in an SNS structure in 
which the quasiparticle current is carried by MAR. An- 
ticipating the experimental shot noise results presented 
in this paper we demonstrate that in a system in which 
multiple Andreev reflections occurs the quasiparticle shot 
noise at V < 2A/e has a maximum value given by 
Si — (1 + |0-)2e/ because current is effectively carried 
by multiply charged quanta. The maximum shot noise 
is obtained if the transmission probability of the system 
approaches zercn. 

The structure under study con- 
sists of a NbN/MgO/NbN SIS tunnel junction with small 
defects in the 1 nm thick MgO barrier acting as paral- 
lel SNS point contacts. Commonly used Nb/A103;/Nb 
tunnel junctions are known to exhibit very low subgap 
"leakage" current because the barrier is formed by ther- 
mal oxidation of a thin Al layer which wets the Nb sur- 
face. For NbN junctions such a thermal oxidation process 
is not yet established. Currently the best results axe ob- 
tained by direct deposition of the barrier materials, and 
defects occur due to "missing atoms" since the barrier is 
comprised of meteiy 1-2 atomic layers in high current 
density junctionalJ. In Fig. |l| the nearly perfect tunnel- 
ing /, V characteristic of a Nb junction is compared with 
that of a typical NbN junction. Both are measured at 
4.2 K. The NbN junction area is 0.8 /xm^, the normal re- 
sistance Rn is 25 J7, the gap voltage 2A/e is 4.8 mV. The 
Nb junction has an area of 1 /.tm^, a normal resistance 
Rn of 56 and a gap voltage of 2.8 mV. Clearly the 
subgap current of the NbN junction is much larger than 
the current in the Nb junction. In addition the subhar- 
monic gap structure in the NbN /, V and dV/dl curves 
indicates thatj-the subgap current is predominantly car- 
ried by MAR.eI. Since MAR is evidence for the presence 
of higher order processesllil the observation of MAR in- 
dicates conducting paths with transmissivities close to 1, 
which we further call pinholes. 
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FIG. 1. Normalized I,V and Rdyn = dV/dl curves of a 
Nb and NbN junction. Tiie subgap structure is clearly visible 
up to the fourth harmonic in the NbN Rdyn curve. 

A remarkable feature of a pinhole system is that vary- 
ing the barrier thickness changes only the number j-of 
pinholes without affecting the pinhole conductance^. 
Therefore the pinholes can be considered as reproducible 
barrier defects with identical transmission probabilities. 
From the relative height of the current steps at the sub- 
harmonic gap structure in Fig. |l| aj^inhole transmis- 
sion probability T of 0.17 is derivedld. The supercur- 
rent (Fraunhofer) dependence on magnetic field is nearly 
ideal, indicating a large number of pinholes with a ho- 
mogeneous distribution over the barrier. All measured 
high current density NbN/MgO/NbN junctions exhibit 
similar behavior. 

The noise measurement setup is schematically depicted 
in Fig. |. 
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FIG. 2. Schematic layout of the current-voltage and noise 
measurement system. All leads going into the dewar are 
low-pass filtered except for the 1.5 GHz wiring. The Al leads 
suffice as heat sink. 

All DC leads are filtered, 30 SI manganin wires carry 
currents between room temperature parts and those at 
cryogenic temperatures. The isolator is used to dissipate 
power reflected off the amplifier into a 50 load in or- 



der to avoid gain variations with varying junction output 
impedance. A band pass filter transmits an 85 MHz band 
centered at 1.5 GHz. Although a spectrum of the signal 
was not taken, we assume the measured noise to be white 
since 1.5 GHz is high enough to ignore 1// noise of the 
device. 

The total gain of the amplification section is measured 
by mounting the high quality Nb SIS junction shown in 
Fig. [| in the measurement setup. Since the Nb junction 
exhibits single particle tunneling, the shot noise gener- 
ated by the current in the juiTjetion can be accurately 
calculated from the /, V curvaia. Corrections due to 
the measurement frequency of 1.5 GHz are negligible for 
voltages V hu K, 0.05mV. The effect of finite tem- 
perature below 0.5 mV is taken into account by using 
Si = 2el coih-{eV /2k bT) to calculate the noise coming 
from the junction. Comparing this noise with the noise 
power after filtering and amplification gives the total gain 
as well as the noise added by the isolator, cables and am- 
plifier. The output power is given bjil3 



.1 



Pout — G, 

+ keil ~ G^so)T) + Pa,np 



ampB { — SlRdyni^ T )Giso 



(1) 



in which Gamp is the amplifier gain, B is the bandwidth, 
/, V and Rdyn are the current, voltage and the differen- 
tial resistance dV/dl of the junction, F is the reflection 
coefficient \Rdyn - Ramp\/{Rdyn + Ramp), Giso IS the iso- 
lator gain, and Pamp stands for the output noise of the 
amplifier chain. From the measured and calculated noise 
power curves the gain and noise of tiie amplifier and iso- 
lator can be accurately determinedli3 using Eq. 0. The 
loss of the isolator plus cables is 0.25 dB. The total am- 
plification including cable losses and reflections is 80 dB. 
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FIG. 3. Contributions of the Nb junction, isolator, and 
amplifiers to the total output noise. The measured total noise 
power is accurately given by Eq. 1. The deviation above the 
gap voltage may be connected to a proximity effect layer near 
the tunnel barrier causing doubled shot noise. 
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The shot noise Si generated in the NbN junction is 
obtained by performing a measurement of Pout, I and 
dl/dV as function of voltage similar to that conducted 
for the Nb junction. Since the loss and gain of each 
component is known, the shot noise Sj can be derived 
from the output noise Pout using Eq. |^. The result is 
plotted in Fig. ^ together with the Poisson shot noise 
2el calculated using the measured /, V curve. The shot 
noise is expressed as a current by dividing the measured 
Si by 2e to allow easy comparison with the shot noise 
calculated form the /, V curve. If the junction behaves 
as an ideal shot noise source as the Nb junction does, 
a comparison of Si/2e with the dc current / will find 
them identical at voltages above 0.5 mV. Clearly, below 
the gap voltage the measured noise is much larger than 
Poisson noise. 
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FIG. 4. Measured and calculated shot noise generated 
by the NbN junction, all normalized to current by division 
by 2e. The square dots give the measured noise, the dashed 
line is the calculated Poisson shot noise curve. The calcula- 
tion of the Si curve labeled (MAR) is explained in the text. 
Above the gap voltage the shot noise follows Poisson noise 
since current transport through the barrier dominates. The 
error bars, given at specific points to avoid crowdedness of 
the figure, refiect the uncertainty in the noise contributions 
of the measurement equipment. Errors at intermediate volt- 
age values can be roughly linearly interpolated. Offsets in 
current and noise are negligible since at negative voltage bias 
the same results are obtained. 

We explain this excess noise by taking into account 
multiple Andreev reflection processes in the pinholes. For 
an NS structure with a disordered region of length L 
which is much smaller than the mean free path, I in this 
region, the shot noise at ^ <C A/e is given by:l 



Pns(T) = ^2 - r)2 ' 



(2) 



in which T is the transparency of the barrier between 
the normal- and superconducting region. For T ^ 1 
Pns — 2PpoiB3on indicating uncorrelated current of parti- 
cles with charge 2e due to Andreev reflection. Realizing 
that via multiple Andreev reflection much larger charge 



quanta can be transferred!, the shot noise generated in 
pinholes by MAR is derived in the following way. The 
pinhole is assumed to consist of two SN and NS structures 
in series. If the measured transmission probability of 0.17 
is inserted into Eq. ^ the shot noise suppression is a mere 
1 %, below the measurement accuracy. Therefore if the 
effectively transferred charge q is known, the shot noise 
is simply found from Si = 2ql. We calculate the charge 
q{V) in the limit of unity pinhole transmission probabil- 
ity using the trajectory method employed in the original 
paper on MARQ. This enables separate calculation of 
the currents Im carried by m Andreev reflections. The 
/,„ values form weight factors used to calculate the av- 
erage transferred charge carried by what might be called 
Andreev clusters. 

The total current flowing from the left superconduc- 
tor (at voltage V) to the right superconductor (at zero 
voltage) is given by the difference oiD 



P 



1 

eRo 



dEfo{E - eV)[l - A{E - eV)]- 



[I + A{E) + A{E)A{E + eV) 
and 



(3) 



dEfo{E)[l-AiE)]- 



[1 + A{E - eV) + A{E - eV)A{E - 2eV) 



(4) 



where Rq is the total normal state resistance of the pin- 
holes, fo{E) is the Fermi distribution and A{E) is the en- 
ergy dependent Andreev reflection probability. The first 
term [1 — A] denotes the fraction of available electrons in 
the superconductor which is transmitted into the normal 
region. The last term [1 + • • •] gives the transferred charge 
multiple which is le if no Andreev reflections take place, 
2e for one Andreev reflection, and so on. The currents Im 
(to = 1,2,...) carried by TO-electron processes are calcu- 
lated by splitting Eqs. ^ and ^ into the TO-electron parts. 
For example I2 is given by 



2 

— / dE{ ME-eV)[l-A{E 
eRo J -00 



,V)] 



A{E)[1 - A{E + eV)]- 



foiE)[l - AiE)]AiE - eV)[l - A{E ~ 2eV)\ ) (5) 

Since the current Im is carried only by carriers with an 
7Ti-fold charge, the magnitude of this current divided by 
the total current gives the relative contribution to the av- 
erage transferred charge. Therefore the average charge q 
of an Andreev cluster for a given voltage V is given by: 



J2m=im- I^njV) 



(6) 
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The resulting charge- voltage curve is shown in Fig. || and 
is used to calculate the noise Sj = 2q{V)I in Fig. ^. 
The measured charge values in Fig. ^ are obtained by 
dividing the measured shotnoise by 2el. The similarly 
derived unity charge of the Nb junction is shown for 
comparison. The effect of finite temperatures on the 
shot noise correction factor coth{eV /2k bT) at voltages 
V < 2A/e is negligible since coth{q{V)V /2kBT) w 

C0th((l + |#).5fy)«C0th(^)«l. 



\\ ' eV ' ZksJ- ' \/crJ / 

Apart from a superposed subharmonic gap structure 
the effective NbN charge is proportional to 1/V. This 
can be appreciated by performing the charge calculation 
ignoring Andreev reflection for energies \E\ > A. Higher 
order terms vanish and the average charge is equal to the 
charge quantum transferred; q„ = ne at voltages Vn — 
-2A. for n = 2,3, . . . giving an analytical approximation 
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FIG. 5. Effective charge as function of voltage. The square 
dots give the charge derived from the measurements. The 
calculated charge values are used to calculate the shot noise 
in Fig. ^. For comparison the effective charge of the Nb cali- 
bration junction is shown by circles. 

In conclusion, for the first time shot noise much larger 
than Poisson noise is observed in an SIS structure. A 
NbN SIS tunnel junction with pinholes is employed to 
obtain current transport dominated by multiple Andreev 
reflections (MAR). Due to the occurrence of MAR large 
charge quanta are transferred between the electrodes, 
causing a significantly enhanced shot noise. A simple 
model is developed to calculate the effective charge from 
which the shot noise is obtained. The model answecs. 
the question about the origin o£_excess noise in NbNO 
and very high current density NbL3 SIS junctions in THz 
applications. 

Note added in proof. After submission j-of this 
manuscript we became aware of theoretical workEj which 
analyzes MAR enhanced shot noise, in particular the de- 
pendence of Si on temperature and transmission. In our 
semi-empirical theory we circumvent the problem of cal- 



culating Si by calculating q{V) (ignoring coherence ef- 
fects) and using the measured /. 
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